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We study the charge and spin transport under high electric field (up to several kV/cm) on the 
surface of topological insulator Bi2Se3, where the electron-surface optical phonon scattering dom- 
inates except at very low temperature. Due to the spin mixing of conduction and valence bands, 
the electric field not only accelerates electrons in each band, but also leads to inter-band precession. 
In the presence of the electric field, electrons can transfer from the valence band to the conduction 
one via the inter-band precession and inter-band electron-phonon scattering. The electron density 
in each band varies with the electric field linearly when the electric field is strong. Due to the 
spin-momentum locking, a transverse spin polarization, with the magnitude proportional to the mo- 
mentum scattering time, is induced by the electric field. The induced spin polarization depends on 
the electric field linearly when the latter is small. Moreover, its magnitude is inversely proportional 
to the temperature and is insensitive to the electron density at high temperature. Our investigation 
also reveals that due to the large relative static dielectric constant, the Coulomb scattering is too 
weak to establish a drifted Fermi distribution with a unified hot-electron temperature in the steady 
state under the electric field. After turning off the electric field in the steady state, the hot carriers 
cool down in a time scale of energy relaxation which is very long (of the order of 100-1000 ps) while 
the spin polarization relaxes in a time scale of momentum scattering which is quite short (of the 
order of 0.01-0.1 ps). 

PACS numbers: 73.50.Fq, 75.70.Tj, 72.25.Rb, 72.25.-b, 71.10.-w 



I. INTRODUCTION 

Topological insulators are a new class of materials that 
attract much interest recently^— They are band in- 
sulating in the bulk but conducting along the surface 
due to the gapless surface state. In three-dimensional 
strong topological insulators, the surface state involves 
an odd number of massless Dirac cones in the sur- 
face Brillouin zone and is protected by time-reversal 
invar iance i&l i 1 2 ' 1 7 ' 1 8 For the simplest case with a single 
Dirac cone, such as on the surface of Bi 2 Se 3 , Sb 2 Te3 and 
Bi2Te3f2~— the surface state near the Dirac cone can be 
described by the Rashba spin-orbit coupling^ 9 - and hence 
exhibits hclicity^ ' 20 ! 21 The metallic and helical surface 
state has been proposed to have much potential for the 
application in spintronics^i 7 - 

The spin helicity of surface state in topological insu- 
lators has been experimentally measured by spin-angle 
resolved photoemission spectroscopy^ 2 - However, it has 
not been clearly evidenced by the electrical transport ex- 
periment. This might be caused by the difficulty in seper- 
ating the bulk and surface conduction as both of them are 
usually involved simultaneously, and also the influence of 
the stray field from the ferromagnetic electrode when the 
spins are injected via the ferromagnetic contact. To over- 
come these circumvents and confirm the spin-momentum 
locking in the topological surface state, a scheme of trans- 
port experiment by injecting spin polarized electrons into 
the topological surface state via silicon has been proposed 
very recently— In spite of the experimental difficulties, 
the understanding of charge and spin transport on the 
surface of topological insulators is necessary. In fact, the- 



oretically, this issue has been preliminarily studied, with 
electric field being small (~0.1 kV/cm ; 12 ' 24 under which 
the electrons are near equilibrium), the Fermi level lo- 
cated deep enough in the conduction band and the fully 
occupied valence band being irrelevant j 1 2 ' 24 ~ ^£ Moreover, 
even in the microscopic study based on kinetic equations 
by Culcer et al, only the electron-impurity scattering 
is considered^ The electron-phonon scattering, which 
can be very important j21 and also the electron-electron 
Coulomb scattering are not incorporated. It is revealed 
that in the presence of driving by the electric field or dif- 
fusion by the density gradient, a transverse spin polariza- 
tion is induced due to the spin-momentum locking; 12 ' 24 ' 25 
with the magnitude proportional to the electric field and 
momentum scattering timei 12 ' 24 Besides, due to the spin- 
momentum locking again, the spin polarization relaxes in 
a time scale of momentum scattering in the absence of 
an electric field i 25 ' 26 

The transport on the surface of topological insulators 
under large electric fields (~1 kV/cm), which can drive 
the carriers far away from the equilibrium. ) 28 ' 29 has been 
rarely investigated so far. Here we perform this study on 
the surface of Bi2Se3 by means of the kinetic spin Bloch 
equation (KSBE) approach ; 30 ' 31 with the electric field 
up to several kV/cm. Bi 2 Se3 is expected to have a 300 
meV direct band gap at the zone center^ and is usually 
n-type due to the charged Se vacancies. However, the 
electron density can be adjusted by counterdoping with 
C a 3 ' 32 ' 33 or partially substituting Bi with Sb to reduce Se 
vacancies.—^ Particularly, employing Cd doping in com- 
bination with a Se-rich growth condition, even a p-type 
Bi2Se3 can be obtained^ For the thin- layer structure 
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with a thickness ~10 nm, the electron density is also 
controllable by the gate voltagej 14 ' 15 Due to the large 
relative static dielectric constant , 14 i 24 i 36 i 37 the dominant 
scattering on the surface is the electron-surface optical 
phonon scattering. 27 

In this work, we take into account both the n-doped 
degenerate and intrinsic nondegenerate cases before turn- 
ing on the electric field, by including both the conduc- 
tion and valence bands. Our study reveals that due to 
the joint effects of the driving of the electric field and 
the inter-band precession as well as inter-band electron- 
phonon scattering, electrons can be transferred from the 
valence band to the conduction one. This effect is more 
pronounced when the scattering is weak and the electric 
field is strong. Moreover, the variation in electron den- 
sity for each band is linear in the electric field when the 
latter is strong for both the n-doped and intrinsic cases. 
The induced spin polarization is linear in the electric field 
when the latter is small and deviates from the linear re- 
lation when the latter becomes large enough. Besides, at 
high temperature, the induced spin polarization is sen- 
sitive to the temperature but insensitive to the electron 
density. We also find that the electron-electron Coulomb 
scattering is too weak to establish a drifted Fermi distri- 
bution with a unified hot-electron temperature under the 
electric field. After turning off the electric field, the hot 
carriers cool down in a time scale of 100-1000 ps while 
the spin polarization relaxes in a time scale of momentum 
scattering, which is of the order of 0.01-0.1 ps. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model and KSBEs. In Sec. Ill, we analyti- 
cally solve the KSBEs under a small electric field in the 
presence of electron-impurity and electron-phonon scat- 
terings, with the latter treated in the elastic scattering 
approximation. In Sec. IV we present the numerical re- 
sults from full calculation based on the KSBEs. We sum- 
marize in Sec. V. 



II. HAMILTONIAN AND KSBES 

We set the z-axis along the (001) direction of Bi2Se3. 
In the collinear spin space spanned by the eigenstates of 
a z , {| t): I I)}, the single-electron Hamiltonian describing 
the low-energy (001) surface state around the T point has 
the Rashba spin-orbit coupling 1 ^ forr" 2 ' 2 ° i21 



Hq = foj(kxz)'<T, 



(1) 



where Vf ~ 5 x 10 5 m/s is the Fermi velocity [vf is exper- 
imentally measured to vary from 2 x 10 5 m/s (Ref. 
to 6 x 10 5 m/s (Ref. @)], k is the two-dimensional elec- 
tron momentum and er is the vector composed of the 
Pauli matrices for spin. The eigenstates of Ho, |k±) = 
-^(±ie~ l6k \ f) + I I)) with # k representing the polar 
angle of k, correspond to the conduction and valence 
bands with linear dispersion £k± = ihvfk. The inter- 
action Hamiltonian H] includes the electron-impurity, 24 



electron-surface optical phonon— and electron-electron 
Coulomb interactions. 

The KSBE o 30 i 31 are applied to study the charge and 
spin dynamics on the surface of Bi2Se3 in the weak scat- 
tering regime where Vf{kr k ) 3> 1 (r^ is the momentum 
scattering time). In the helix spin space constructed by 
{|k+), |k-)},2L the KSBEs reac^i 



d t pk(t) +i[vfka z ,pu{t)} - {eE/K)d kx Pk{t) 
- (eE/H)lUtd k „U k ,p*(t)]+d t p k (t)\. = 0. 



(2) 



Here pu(t) is the 2x2 density matrix for electrons 
with momentum k in the helix spin space. The diag- 
onal elements of Pk(i), P\t++/--(t) = fk+/-(t)-> stand 
for the distributions of electrons in the conduction and 
valence bands respectively, and the off-diagonal ones, 
PkH (t) = p k i-C0i characterize the inter-band coher- 
ence. The second term on the left-hand side of Eq. (|2J| is 
the coherent term. The third and fourth terms are con- 
tributed by the electric field along the x-axis, E = Ex. 
The third term accelerates electrons in each band, while 
the fourth term, with 



(3) 



and hence Uld kj; U k = ts ' 2 n fc 8k (1 - a x ), leads to the inter- 
band precession. This electric-field-induced inter-band 
precession originates from the spin mixing in the con- 
duction and valence bands. In fact, this effect also ex- 
ists in graphene where the pseudo-spins are mixed in the 
conduction and valence bands, as revealed by Balev et 
aJ.j2& However, in their study with a low electric field 
~0.01 kV/cm, the inter-band coherence and hence nat- 
urally the precession term are neglected in their kinetic 
equations. 38 Here we retain this term in the KSBEs and 
will show it influences the transport properties markedly. 
The last term on the left-hand side of Eq. @ is the scat- 
tering term, explicitly given in Appendix \X\ Due to the 
large relative static dielectric constant kq of Bi2Se3 (we 
typically set kq — 100 following Refs. [lj and while 
kq w 50 ~ 200 as given in the literatur o 36 ' 37 ) and the 
low energy of surface optical phonons (h&jQ — 7.4 meV), 
the electron-electron and electron-impurity scatterings 
are less important than the electron-phonon scattering 
in a large temperature regioni 37 i 39 ' 40 

By solving the KSBEs in the presence of electric field, 
we can investigate the charge and spin transport proper- 
ties on the surface of topological insulator. The electron 
and hole densities in the two bands are n e — ^ k /k+ and 
Uh = 2k(l — /k-)j respectively The spin polarization 
in the collinear spin space, induced by the electric field, 
reads 



s = Es k = ^E Tr M 



(4) 



where pk — U^p^U^ is the electron density matrix in 
the collinear spin space. Due to the Rashba spin-orbit 
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coupling [Eq. (JTJ], the charge current density along the 
x-axis under the electric field reads 



evf ^2 Tr[p k a y ] = 2ev f S y /h, 



(5) 



exactly proportional to the induced spin polarization. 
Starting from the steady state established by the elec- 
tric field, we can study the cooling of carriers and the re- 
laxation of spin polarization by solving the KSBEs in the 
absence of the electric field. In the following, we carry out 
these studies first analytically under a low electric field 
and then numerically in the large electric field regime. 



III. ANALYTICAL STUDY WITH ELASTIC 
SCATTERING AT LOW ELECTRIC FIELD 

Under a small electric field (eE(r k ) <C h{k)), the KS- 
BEs can be solved analytically with only the elastic scat- 
tering, e.g., the electron-impurity scattering. Due to the 
low energy of surface optical phonons, we can also incor- 
porate the electron-phonon scattering by elastic scatter- 
ing approximation. Therefore, with both the electron- 
impurity and electron-phonon scatterings, the scattering 
term in Eq. @ can be simplified from Eq. (jAll) as 



mentum scattering rates 
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(9) 
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/ d0HAj| 2 + (27Vo + l)| 5q | 
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x sin0 sinW. 



(10) 



To obtain Eq. ([8]), we have replaced Pk{t) in the terms 
proportional to E in Eq. ([2]) by the initial equilibrium 
state jOk(0) = diag{/k+(0), /k-(0)} when E is small, 
and also assumed |-D q | 2 to depend only on the magni- 
tude but not the direction of q (q = 2fcsin|) by treat- 
ing the screening in the long-wavelength and static limit 
for analytical feasibility (refer to Appendix [SJ. Here 
/k±(0) = l/{exp[/3(ek± — m)] + 1} i s the Fermi distri- 
bution, where ft — l/(fcsT) and /i is the chemical po- 
tential. In both the degenerate and nondegenerate lim- 



its, \D q 



approximately (refer to Appendix [XJ 



and hence the momentum scattering rate limited by the 
electron-impurity scattering is cx fc _1 . However, for the 
electron-phonon scattering, with |g q | 2 sa \g \ 2 = |<? q= o| 2 , 
the momentum scattering rate limited by it is cx k. 



n-2n 

\ d0k'MAc-kf + (2iVo + l) 
Jo 



87rfi. 2 u/ 

x | ffk -k< | 2 ] {2[1 + cos(0 k - 0k')] [pk(t) ~ pv(t)] 
+ isin(0 k -0 k O[pk(t)-/v(*)>^]} • ( 6 ) 

J k' — k 

Here |-Dk-k'| 2 and |<?k-k'| 2 are the electron-impurity (re- 
fer to Appendix [XJ and electron-phonon scattering ma- 
trix elements respectively, rii is the impurity density, and 
N a is the number of optical surface phonons. Explicitly, 



ISql 



(Ai + \ 2 q) 2 h/(2Amu JQ ), 



(7) 



where M is the ion mass, 21 is the primitive cell area 
[ft 2 /(M2l) = 4 x 1CT 3 meV], Ai « 5 eV-nm and A 2 « 
1.6 eV-nm 2 i22 In our study, the largest q involved is of 
the order of 0.1-1 nm _1 and therefore the momentum 
dependence of |g q | 2 is weak in the presence of the large 
constant Ai. It is noted that in Eq. ([6]) the contribution 
with 1 0k — 0k' | = 7T is zero, indicating the absence of 
backscattering in the same band. 

By performing the Fourier transformation on Eq. @ 
with respect to the polar angle 0k and defining 8 k (t) = 

r 2ir 



h Jo d0^e- Ud -TT[y k (t)a}, one has 



d t S l k (t) + 2v f kS{(t) x z + S l k (t)/4 k - iS l k (t) x x/r, 
- e£^, ± iAfeZ/4 - ieElSi^Alk^y/A = 0. 



I 

2k 

(8) 



Here A° k 



/k+(0) - /k-(0). 



9fcA? and the mo- 



A. Spin polarization induced by the electric field 



Retaining the lowest three orders of S' fc in Eq. ©, i.e., 
those with I = and ±1, one can obtain the steady-state 
solution as 



Sg(+oo) = (0,0,fiAg/2) 
1±u , \ TieE 



0+°°) 



StH+oc) 



ky 



s£ x (+°o) 



8v f k 

±ieE 
16(vfk) 2 Tk 
eEr k ^ 



(Alk-' + Al), 



Alk- 1 



Afc), 



A^ 



4(vfkr k ) 



(11) 
(12) 

(13) 

(14) 



where r k = 
with Vf(kr k ) 



Sf (+oo) 



rl k = r^ k - In the weak scattering limit 
S> 1, S^ 1 (+cxj) can be approximated as 
eEr k A k z/ 4. The steady-state solution 



^kz(^~°°) ~ ^fez(O) means that the electron density in 
each band keeps unchanged under the small electric field 
in the framework of this analytical study. 

Now we focus on the steady-state spin polarization in 
the collinear spin space. According to Eq. (U), the spin 
polarization reads S = dkkS k with 

S° = (-ImSl + ReSl y , -ImS^ - ReS^, -S° kx ). (15) 

With the solution of S"' 1 (+oo) obtained above, we have 



S°(+oc) « -eEr k A\y/A. 



(16) 
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Therefore, the total spin polarization S(+oo) ~ 
e£;(T fe )ln(2 + + e~^)y /(8nphv f ) where (r fe ) = 

f + °° dkkAl.Tk/ f °° dkk/S.\. For the intrinsic nondegen- 
erate case with fi = 0, S(+oo) w eE(rk) ln4y/(8n/3Hvf), 
and for the n-type degenerate one, S(+oo) rs 
eETfc, kfy /(87t) where kf is the Fermi momentum. The 
spin polarization obtained for the degenerate case is con- 
sistent with that presented in Ref. [l2| where the valence 
band is fully occupied. These results indicate that under 
the small electric field, a transverse spin polarization is 
induced due to the Rashba spin-orbit coupling [Eq. (JTJ], 
with the magnitude proportional to the electric field and 
momentum scattering time. 

We further study the momentum scattering time in 
the intrinsic nondegenerate and n-type degenerate cases 
respectively, with the help of Eqs. (|9j)- (|10[) . For the 
intrinsic case, the mean momentum scattering times 
limited by the electron-impurity and electron-phonon 
scatterings are (r^ 1 ) rs 8(3hn e / (nr^m In 4) and (r^ p ) w 
4/3ft 3 i> 2 -/[(27V + l)|g | 2 ln4], respectively. For the n- 

type case, at the Fermi level t" kf /[wr^VfniI(r a )] 

and r£ P rs 4h 2 v f /[(2N + l)\go\ 2 k f ]. Here I(r s ) = 

Jq 27 "^ sin 2 0/(2 sin | + r s ) 2 and r s = e 2 / (47re K hvf), 
with eo being the static dielectric constant [when kq = 
100, one has r s rs 0.044 and I(r s ) rs 2.6]. To calculate 
the momentum scattering time limited by the electron- 
phonon scattering, we have neglected the weak momen- 
tum dependence of |<? q | 2 by approximating it as |<?o| 2 . 

Consequently, the induced spin polarization lim- 
ited by the electron- impurity scattering is S cl (+oo) rs 
eEn e y / (ir 2 r 2 v frii) for the intrinsic nondegenerate case 
and S CI (+oo) rs eEn e y /[2TTVfr 2 I(r s )rii] for the n- 
type degenerate one. However, the one limited 
by the electron-phonon scattering is S ep (-(-oo) rs 
eEh 2 vfy /[2n\go\ 2 (2N + 1)], in spite of the electron den- 
sity. Considering both scatterings, the induced spin po- 
larization along the y-axis reads 



S„(+oo) « [l/S^+oo) + l/S p (+oo)] 



E ■ 10 12 cm/kV 
ai(27V + 1) + a 2 ni/n, 



-h ■ cm 



(17) 



with ax rs 13 and a-i rs 0.6 (1) for the intrinsic nonde- 
generate (n-type degenerate) case. It is interesting to see 
that, when the electron-phonon scattering dominates, the 
induced spin polarization is only sensitive to the temper- 
ature. Especially, at high temperature T » hcuo/kB ~ 
90 K, one has S y (+oo) oc T _1 . Due to the large relative 
static dielectric constant, only when the impurity density 
is high enough (e.g., with n, > 10n e ) and the tempera- 
ture is low, the effect of electron-impurity scattering can 
be comparable to that of the electron-phonon scattering 
and then leads to the impurity /electron density depen- 
dence of spin polarization. 



B. Spin relaxation 

We then start from the steady state obtained previ- 
ously to study the spin relaxation with the electric field 
turned off. From the initial state S°' 1 (0) = S^(0)z, one 
can solve the temporal evolution of S fe ' (t) by Eq. ([5]) 
after the terms proportional to E have been removed. 
Then according to Eq. ([15]). one obtains S£(i) = -§L(t)y 
with 



S° ky (t) = SL(0)e 



c k t 



Ckt, 



z k (1 — Cfesin cos — ) 



Tk 



(18) 



Here S° ky (0) = -5^,(0) and c k = y/(2v f kT k ) 2 - 1. 
In the weak scattering limit with Cfe » 1, S ky (t) rs 

^(0)e" t/rfc , indicating that the spin polarization re- 
laxes in the time scale of momentum scattering. This 
feature is in agreement with that given by Schwab et 
aL— and Burkov and Hawthorn.— 



IV. NUMERICAL RESULTS 

The analytical study in the previous section only ap- 
plies to the small electric field. In order to take into 
account the large electric field, as well as all the scat- 
terings explicitly, we carry out the numerical calculation 
based on the KSBEs. In the calculation, we first apply 
the electric field along the x-axis to study the charge and 
spin transport, and then turn off the electric field after 
reaching the steady state to look into the cooling of the 
hot carriers and the relaxation of the previously induced 
spin polarization. We consider both the intrinsic non- 
degenerate and n-doped degenerate cases, starting from 
the initial equilibrium state with fi = and a given n e (0) 
respectively. In our study n e (0) is chosen to be of the or- 
der of 10 11 cm~ 2 or even smaller, to avoid entering the 
bulk states in the presence of high electric field and also 
ensure the validity of Hamiltonian Ho [Eq. (JTJ] without 
involving the terms square or cubic in moment um j 2 1 1 24 
We also do not consider the intrinsic case under low tem- 
perature as for such case the carrier density is quite low 
and the fluctuation (e.g., the effect of puddles 14 ), which 
is beyond the scope of this work, becomes important. 



A. Redistribution of electrons between two bands 

In semiconductors with a large band gap, the electron 
and hole densities in the conduction and valence bands 
remain unchanged under the static electric field. How- 
ever, here on the surface of topological insulator, due 
to the spin mixing of the two bands, the static electric 
field leads to inter-band precession^ 38 and also, due to 
the zero band gap, the inter-band electron-phonon scat- 
tering is easy to take place. Therefore, electrons can be 
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transferred from the valence band to the conduction one 
under the influence of the electric field, with the differ- 
ence in the densities of electrons and holes in the two 
bands, n e — n/,,, keeping constant due to the particle con- 
servation. This redistribution of electrons between two 
bands under the electric field has also been revealed in a 
similar system, the gapless graphene^ In principle, the 
Auger process of the Coulomb scattering, during which 
one and only one of the two scattered electrons transfers 
between two bands ; 41 ' 42 may also lead to the redistribu- 
tion of electrons between two bands under the electric 
field. However, this process is actually forbidden when 
the dynamic screening under random phase approxima- 
tion (RPA^ii&ii (refer to Appendix EJ is adopted.- 4 ^ 
The redistribution of electrons in the two bands can not 
be revealed by the analytical study presented in Sec. IIII1 
which fails to incorporate the inter-band precession [due 
to the replacement of pk(t) by Pk(0) in terms propor- 
tional to E] as well as the inter-band scattering (due to 
the elastic scattering approximation). However, this pro- 
cess can be revealed by fully solving the KSBEs, as stated 
in the following. 

We take an n-type degenerate case with T = 50 K 
and n e (0) = 10 10 cm~ 2 (the corresponding Fermi en- 
ergy is £f ~ 12 meV and the Fermi momentum is 
kf « 0.035 nm _1 ) as an example to show the tempo- 
ral evolution of the relative change of electron density in 
the conduction band, n e (t) / n e (0) — 1, under electric field 
E = 0.1 kV/cm. It is seen in Fig. [TJa) that when both 
the inter-band precession and inter-band electron-phonon 
scattering are excluded, n e (t) remains almost unchanged 
(chain curve). However, once the inter-band precession is 
included, electrons can be effectively transferred from the 
lower band to the higher one (dotted curve). The scenario 
is that, the electrons are precessed from the states below 
the Dirac point to the ones above and then driven away 
to higher-energy states by the electric field. The transfer 
of electrons can also be alternatively realized with the 
assistance of phonons (dashed curve). Around the Dirac 
point (with k < ko = u>o/vf w 0.022 ran -1 ), when elec- 
trons in the conduction band are driven away to higher- 
energy states, electrons in the valence band tend to enter 
the conduction band by absorbing phonons. However, the 
above two effects are not definitely superimposed when 
both of them are present. When the inter-band pre- 
cession transfers electrons between two bands effectively, 
i.e., leads to substantial noncquilibrium between the two 
bands (such as the case presented here), electrons in the 
conduction band tend to fall back to the valence band 
by emitting phonons. Therefore, when both the inter- 
band scattering and inter-band precession are included, 
the steady-state electron density in the conduction band 
decreases compared to the case with only the inter-band 
precession (compare the solid and dotted curves). 

As pointed out above, the inter-band precession and 
inter-band scattering open channels for electron transfer 
between the two bands. However, the steady state of 
the two-band system is determined by the balance be- 
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FIG. 1: (Color online) Temporal evolution of n e (t)/n e (0) — 1 
with E = 0.1 kV/cm, T = 50 K and n e (0) = 10 10 cm" 2 . 
During calculation, the electron-impurity scattering is not in- 
cluded while the intra-band electron-phonon scattering (la- 
beled as intra-ep) is always present, (a): The chain curve 
is calculated without the inter-band precession (labeled as 
IP) and inter-band electron-phonon scattering (labeled as 
inter-ep), the dashed (dotted) one with only the inter-band 
electron-phonon scattering (inter-band precession), while the 
solid one with both the inter-band precession and inter- 
band electron-phonon scattering. The double-dotted chain 
curve is calculated with the inter-band precession, inter-band 
electron-phonon scattering and also the electron-electron 
Coulomb scattering (labeled as ee). (b): Temporal evolu- 
tion of n e (t) /n e (0) — 1 with the electron-phonon scattering 
artificially strengthened. The inter-band precession is always 
included. The solid curve is calculated with both the intra- 
and inter-band electron-phonon scatterings strengthened by 
a factor Xc P = 2, while the dotted (double-dotted chain) one 
with only the intra-band (inter-band) electron-phonon scat- 
tering strengthened by a factor 2. The crosses (with scales 
on the top and right-hand side of the frame) are the relative 
variation of n e reached in the steady state, n e (-|-oo)/n e (0) — 1, 
against the modulation factor Xcp ranging from 1 to 4. The 
chain curve (with scales on the top and right-hand side of the 
frame) is the function [n e (+oo)/n e (0) — l]|x cp =i/Xe P - 



tween the rates of energy gain from electric field and 
the energy loss to phonons. The former is determined 
by d t ei(t) = Ej x (t) cx ES y (t) [refer to Eq. ©]. If the 
steady state is not far away from the equilibrium, one ap- 
proximately has <9 t £i(+oo) cx E 2 \g \~ 2 in the presence of 
electron-phonon scattering only. The dominant channel 
of the energy loss is the intra-band electron-phonon scat- 
tering. That is because the inter-band electron-phonon 
scattering is limited to a small finite region k < ko in 
momentum space, and also, with such small momentum, 
the scattering is weak as the rate is cx k. A rough es- 
timation gives that the rate of energy loss due to the 
intra-band electron-phonon scattering for the degener- 
ate case with both low temperature and electron density 
is d t e (t) cx |(7o| 2 [™ e (£) — n e (0)] (refer to Appendix |B|) . 
Therefore in the steady state not far away from the equi- 
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librium, one has 

n e (+oo)-n e (0)cx£ 2 | 5 or 4 , (19) 

required by d t ei(+oo) = d t e (+oo). 

Based on Eq. (fl"9|) . it is found that the strengthening 
of the electron-phonon scattering, especially the intra- 
band part, leads to the decrease in n e (+oo)/n e (0) — 1. In 
Fig. QJb) we numerically verify this by performing sim- 
ilar calculation as in (a) but with the electron-phonon 
scattering artificially strengthened. It is shown by the 
solid curve that when the total electron-phonon scatter- 
ing is strengthened by a factor Xcp = 2 (solid curve), 
?i e (+oo)/n e (0) — 1 substantially decreases compared to 
the genuine case [solid curve in (a)]. When only ei- 
ther the intra- or inter-band electron-phonon scattering 
is strengthened (dotted and double-dotted chain curves, 
respectively), n e (+oo)/n e (0) — 1 also decreases. How- 
ever, the decrease is not obvious when only the inter- 
band electron-phonon scattering is strengthened (double- 
dotted chain curve), indicating that the contribution to 
the energy loss by the inter-band electron-phonon scat- 
tering is indeed relatively weak. We further show the 
steady-state value, n e (+oo)/n e (0) — 1, against x op rang- 
ing from 1 to 4 by the crosses (the scales are on the top 
and right-hand side of the frame) . It is seen that approx- 
imately n e (+oo)/n e (0) — 1 oc Xcp > satisfying Eq. (Til)]) 
[as a guide to the eye, a function proportional to Xcpi 
[n e (+oo)/n e (0) — l]| Xop =i/Xcpi is plotted by the chain 
curve in (b)]. 

At last we briefly address the effect of electron- 
impurity and electron-electron Coulomb scatterings on 
the redistribution of electrons. Although neither of them 
leads to energy loss directly, both of them can limit the 
charge current and hence the energy injection rate of the 
electric field. The effect of impurities on charge current 
is apparent, as also explicitly indicated here by the com- 
bination of Eqs. (|T7f and ([5]). The Coulomb scattering 
is usually deemed to preserve the charge current. This 
is indeed the case in semiconductors with parabolic en- 
ergy spectrum where both momentum and current are 
conserved during the Coulomb scattering. However, here 
with linear energy spectrum, only the momentum but 
not the current is conserved during the Coulomb scat- 
tering. This particular feature also exists in graphene 
with linear dispersion as well4£ In fact, the effect of the 
Coulomb scattering can be conjectured based on our ana- 
lytical study, which indicates that j x oc S y oc tj~ [Eqs. (0 
and (|16p]. With the addition of the Coulomb scatter- 
ing, the momentum scattering time is reduced and 
hence the charge current decreases. Consequently, with 
the inclusion of the electron-impurity and/or electron- 
electron Coulomb scattering, the current and hence the 
energy injection rate of electric field decreases, result- 
ing in less obvious redistribution of electrons. Neverthe- 
less, due to the large relative static dielectric constant, 
the effect of the electron-impurity and electron-electron 
Coulomb scatterings is expected to be weak except when 
the temperature is low. In Fig. HJa), we add the double- 



dotted chain curve calculated with the Coulomb scatter- 
ing included. In the absence of impurities and under low 
temperature, the contribution of the Coulomb scattering 
is visible, with which the variation in n e decreases com- 
pared to the Coulomb scattering-free case (solid curve 
there). 

B. Charge and spin transport 

In the following, we systematically investigate the 
charge and spin transport under the electric field, first 
in the low electric field regime and then the large one. 

1. Low electric field regime 

We first focus on the low electric field regime with E < 
0.1 kV/cm and compare the numerical results with the 
analytical study in Sec. IIIII We consider three cases, (I) 
n-doped case with T = 50 K and n e (0) = I0 10 cm" 2 , (II) 
n-doped case with T = 300 K and ra e (0) = 5 x I0 11 cm" 2 
and (III) intrinsic case with T = 300 K and /i = [cor- 
respondingly n e (0) « 0.8 x 10 11 cm" 2 ]. Corresponding 
to these three cases, in Figs. [2Ia)-(c), we show the de- 
pendence of n e (+oo)/n e (0) — 1 on the electric field with 
different impurity densities (note that the y-axes are in 
different scales in the three figures). For case (I) with 
low temperature T — 50 K, n e increases obviously with 
the electric field, due to the small initial electron den- 
sity in the conduction band n e (0) as well as the weak 
electron-phonon scattering. Moreover, n e increases with 
E faster with a smaller impurity density as the electron- 
impurity scattering weakens. Besides, it is shown that 
when E is small (E < 0.05 kV/cm), n e (+oo) /n e (0) - I 
roughly exhibits a square dependence on E, consistent 
with Eq. (Ti"9"l) . For case (II), which is highly degenerate 
and with strong electron-phonon scattering due to both 
high temperature and electron density, the electron den- 
sity remains almost unchanged in the low electric field 
regime under investigation (the relative variation is of 
the order of 10 -3 ). Nevertheless, for the intrinsic case 
(III), it becomes relatively easier for electrons to transfer 
due to the less occupancy of electrons in the conduction 
band, when compared to case (II). Finally, Figs. [U[b) 
and (c) show that the effect of impurities is negligible at 
T = 300 K. 

The induced spin polarizations in the steady state 
S y {+oo) against the electric field E for the three cases (I)- 
(III) are plotted in Figs.[UJd)-(f), respectively. It is shown 
that S y (+oo) increases with E much faster at 50 K than 
those at 300 K, mainly due to the weaker electron-phonon 
scattering. For case (I) with T = 50 K, the electron- 
impurity scattering is important and n e increases with 
E effectively [refer to Fig. [2a)], therefore the rate of the 
increase in S^+oo) with E shows impurity density de- 
pendence and also a slight electric field dependence [re- 
fer to Eq. (|TT)) ]. However, for cases (II) and (III) with 
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E (kV/cm) 

FIG. 2: (Color online) (a)-(c): The relative change of electron 
density in the conduction band, n e (+oo)/n e (0) — 1, against 
the electric field E for cases (I)- (III), respectively, with differ- 
ent impurity densities, no = 10 10 cm -2 , (d)-(f): The induced 
spin polarization S y (+oo) against the electric field E for cases 
(I)- (III) with different impurity densities (the scale is on the 
right-hand side of the frame). 



T = 300 K and hence the strong electron-phonon scat- 
tering, S y (+oo) increases with E linearly, with almost 
the identical rate under different impurity and electron 
densities. This is consistent with the analytical study: 
when the electron-phonon scattering dominates, the spin 
polarization increases with the electric field linearly in a 
rate solely determined by temperature [refer to Eq. (|17p 
and the discussion there] . 

In the following, we compare the induced steady-state 
spin polarization numerically obtained to the analytical 
one at E — 0.01 kV/cm, with which all the three cases 
(I)-(III) are close to the equilibrium. In Fig. |31 the im- 
purity density dependence of the spin polarization from 
the numerical calculation with all the scatterings explic- 
itly included is shown by triangles for the three cases. 
We also plot the numerical results without the Coulomb 
scattering by crosses. It is shown that the effect of the 
Coulomb scattering can only be visible when the impurity 
density approaches zero and the temperature is low. For 
comparison, the analytical calculation based on Eq. (fTT)) 
is plotted by the solid curves in the figure. One finds that 
the analytical formula primarily captures the numerical 
results. We point out that at 300 K where the elastic 
scattering approximation for electron-phonon scattering 
is more reasonable (fcgT 3> Hu)q), the best fit to the nu- 
merical results by Eq. (1171) requires ol\ = 17. Finally, 
Figs. [2Jb) and (c) indicate more clearly that at room 
temperature the electron-impurity scattering is negligi- 
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FIG. 3: (Color online) Electric-field-induced spin polariza- 
tion S y (+oo) in the steady state against impurity density rii 
for cases (I)- (III), in (a)-(c) respectively. The triangles are 
obtained by numerically solving the full KSBEs, while the 
crosses are calculated without the electron-electron Coulomb 
scattering (labeled as ee). The solid curves are from the ana- 
lytical result [Eq. (TIT)) ]. 

ble even when n» reaches 2 x 10 11 cm" 2 . In fact, accord- 
ing to the analytical study in Sec. IIII A[ the momentum 
scattering times are estimated to be t^ 1 — 4.5 ps and 

rP = 2.5 ps for case (I), r|' — 31.8 ps and = 0.07 ps 

for case (II), and (t^ 1 ) = 19.7 ps and (r^ p ) = 0.17 ps for 
case (III) (in the estimation the impurity density is set 
as rii = 10 11 cm -2 ). These values quantitatively support 
the dominance of the electron-phonon scattering at high 
temperature. 



2. Large electric field regime 

After investigating the low field regime and comparing 
with the analytical study, we proceed with the large elec- 
tric field regime (E is upto 7 kV/cm). Both the intrinsic 
nondegenerate and n-type degenerate cases under 150 K 
and 300 K are considered. For the degenerate case, we 
set n e (0) = 5 x 10 11 cm~ 2 . The impurity density is fixed 
at rii — 10 11 cm~ 2 , with which the electron-impurity 
scattering is in fact negligible under the temperature in- 
vestigated. 

We first present the electric field dependence of the 
steady-state electron density in the conduction band, 
n e (+oo), for the intrinsic and n-type cases at different 
temperatures. Figure @] indicates that n e (+oo) increases 
with E almost linearly for both cases, except in the low 
electric field region E < 2 kV/cm for the degenerate sit- 
uation where the Pauli blocking is important. In fact, 
a nearly linear increase of electron density in conduc- 
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tion band with electric field is also observed in intrin- 
sic graphene when E > 0.01 kV/cnn 3 ® It is also noted 
that here in the nearly linear regime, the rate of the in- 
crease in rt e (+oo) with E is determined by temperature 
and is insensitive to the electron density. It is believed 
that this linear relation as well as the temperature de- 
pendent increasing rate are attributed to the dominant 
electron-phonon scattering with almost constant scatter- 
ing matrix element. As a comparison, we recalculate the 
intrinsic case at 300 K by artificially setting Ai = and 
A2 = 20 eV-nm 2 in the electron-phonon scattering ma- 
trix element [Eq. ([?}]. The result is shown in the figure 
by dots. It is seen that it deviates from the linear rela- 
tion obviously. Finally, we also present the calculation 
without the electron-electron Coulomb scattering for the 
intrinsic case at 150 K by the closed squares. The com- 
parison between the closed and open squares indicates 
that the influence of the Coulomb scattering is marginal 
when the electric field is low and relatively effective when 
the electric field is high. 
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FIG. 4: (Color online) Electric field dependence of the steady- 
state electron density in the conduction band n e (+oo) for the 
intrinsic nondegenerate (solid curves) and n-type degenerate 
(dashed curves) cases under different temperatures. The dots 
are calculated for the intrinsic case at T = 300 K by setting 
Ai = and A2 = 20 eV-nm 2 in the electron-phonon scattering 
matrix element [Eq. (J7J]. The closed squares are calculated 
without the electron-electron Coulomb scattering for the in- 
trinsic case at 150 K. 

We then look into the steady-state spin polarization 
S y (+oo) induced by the electric field. In Fig. [5] the elec- 
tric field dependence of S y (+oo) for the intrinsic and n- 
type cases at different temperatures is plotted. Strik- 
ingly, in this whole large electric field region with the 
electron-phonon scattering being dominant, S y (+oo) is 
solely determined by temperature, as previously revealed 
in the low electric field regime. Moreover, with the same 
electric field E, S v (+oo) at 150 K is about 2 times as 
large as that at 300 K, approximately satisfying the an- 
alytical relation S y (+oo) oc T~ x at high temperature. 
In the figure we also plot the S v (+oo)-E relation given 
by Eq. CE7J) (with modified at = 17) at 150 K (dot- 
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FIG. 5: (Color online) Electric field dependence of induced 
spin polarization in the steady state for the intrinsic (solid 
curves) and n-type degenerate (dashed curves) cases at dif- 
ferent temperatures. The dotted and chain curves are ob- 
tained from Eq. |T7J) with ai = 17 at 150 K and 300 K, 
respectively. The closed squares are calculated without the 
electron-electron Coulomb scattering for the intrinsic case at 
150 K. 



ted curve) and 300 K (chain curve), respectively. It is 
seen that the analytical formula for S y (+oo) can ap- 
proximately apply up to E c — 2 kV/cm (4 kV/cm) at 
150 K (300 K). When E exceeds E c , the electron heat- 
ing becomes important. Therefore, with the occupation 
of larger- momentum states, the electron-phonon scatter- 
ing is strengthened. At the same time, the momentum 
dependence of the electron-phonon scattering matrix el- 
ement [arising from the term A29 in Eq. ([7])] also plays a 
role and further enhances the electron-phonon scattering. 
As a result, S y (+oo) tends to decrease and hence deviates 
from the linear relation against E. Finally, the result ob- 
tained without the electron-electron Coulomb scattering 
for the intrinsic case at 150 K is also plotted by the closed 
squares. Again, the effect of the Coulomb scattering on 
spin polarization is also shown to be marginal, especially 
in the low electric field regime. 

We now turn to study the mobility of the two-band 
system. According to Eq. |(SJ), the steady-state charge 
current under the electric field is immediately obtained as 
j x (+oo) = 2eVfS y (+oo)/h. Therefore, the steady-state 
mobility of the electron-hole system can be determined 
by = jx(+oo)/[eE(n e (+oo) +n h (+oo))]. In Fig.|5]we 
plot \x c against E for the intrinsic and n-type cases at 
different temperatures. It is shown that fi c is around the 
order of 10 3 cm 2 /(V-s), in consistence with the experi- 
mental data^ i 14 ' 33 /j, c decreases with E as electrons and 
holes are heated to occupy large-momentum states where 
strong electron-phonon scattering takes place. Moreover, 
with the increase of temperature from 150 to 300 K, 
/i c decreases as well. In fact, we roughly have fi c ~ 
(fi® 1 + jE)~ 1 as both n e (+oo) + n^(+oo) and ja;(+oo) 
are close to linear functions of E. Here 7 is determined 
by the increasing rates of n e (+oo) and j x (+oo) with E, 
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both of which are only sensitive to temperature when E 
is large (as indicated by Figs. 0] and [5]). Therefore, in the 
large electric field regime, [i c ~ 7~ 1 _E _1 , with 7 deter- 
mined by the temperature. This feature is manifested in 
Fig. E] when E > 3 kV/cm. 
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FIG. 6: (Color online) Electric field dependence of mobility 
for the electron-hole system in the intrinsic (solid curves) and 
n-type (dashed curves) cases at different temperatures. 



C. Effect of the Coulomb scattering 

Now we turn to see the effect of the Coulomb scat- 
tering on the steady state established in the presence of 
electric field. It has been revealed in the previous sec- 
tion that due to the large relative static dielectric con- 
stant, the influence of the Coulomb scattering is marginal 
in the presence of electron-phonon and electron-impurity 
scatterings. With the weak Coulomb scattering, elec- 
trons in two bands fail to reach the drifted Fermi dis- 
tribution with a unified hot-electron temperature in the 
steady state. To reveal this, we take the intrinsic case 
with T = 300 K and E = 3 kV/cm as an example to 
plot the dependence of function = log(l//k+ — 1) 

on k along the x-axis in Fig. [7J with different Coulomb 
scattering strengths adjusted artificially^ The compar- 
ison among these sets of data indicates that only when 
the Coulomb scattering is strong enough (e.g., with the 
Coulomb scattering term rescaled by a factor Xcc = 50), 
in the steady state the drifted Fermi distribution /k± = 
I /{exp[(?ii>/|k — — /i±)/(fcflr e )] + 1} is reached. Here 
T e is the unified hot-electron temperature and is t ne 
shift of the momentum center limited by the scattering. 
Both can be obtained from the slope of the wings of the 
"V" shape and the position of the valley, respectively 
(refer to the solid curve as a guide to the eye). n± are 
chemical potentials for the two bands. For the intrinsic 
case as presented here, /i+ = — due to the symme- 
try between the two bands. In the inset of the figure we 
plot both (closed circles) and ^k- (open circles) 

for the genuine case (with x oc = 1) in a small momen- 
tum scale. Although the drifted Fermi distribution is not 
established, the symmetry between the conduction and 



valence bands is clearly indicated by the inset. 
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FIG. 7: (Color online) Function = log(l// k+ — I) against 
momentum k along the :r-axis in the steady state, evolving 
from the intrinsic case with T = 300 K and E — 3 kV/cm. 
The dots are obtained from the calculation for the genuine 
case, while the open squares, triangles and closed squares are 
from the calculation with the Coulomb scattering term mul- 
tiplied by Xcc = 0, 10 and 50, respectively. The solid line 
is plotted as a guide to the eye by fitting the closed squares 
with hot-electron temperature T e w 488 K. Inset: Function 
jFkT? = logU/Zkn? — 1) [closed (open) circles are for electrons 
in conduction (valence) band with r\ = + ( — )] against mo- 
mentum k along the a;-axis in a small momentum scale in the 
steady state for the genuine case. 



D. Spin relaxation 

Finally, we start from the steady-state spin polariza- 
tion reached previously under the electric field, to study 
the spin relaxation by turning off the electric field. The 
zero point of time, t = 0, is reset to be the moment 
at which the electric field is turned off. Beginning at 
t — 0, the heated electrons cool down to the original 
Fermi distribution in the time scale of energy relaxation, 
which is as long as ^100-f000 ps, mainly due to the 
weak inter-band electron-phonon scattering. However, 
the spin polarization relaxes in the time scale of momen- 
tum scattering which is quite shorti 25 ' 26 According to 
Eq. (fl~8|) . Sk y {t) relaxes approximately in the exponential 
form when Vfkru S> f . The total spin polarization S y (t), 
containing a summation of Sk y (t) over momentum, does 
not definitely relax in a fine exponential form. Our calcu- 
lation shows that S y (t) against t decays first very fast to 
about 0.15^(0) but then slowly. The slowly-decaying tail 
with small magnitude is determined by the low-energy 
states near the Dirac point with long momentum relax- 
ation times. As an example, in the inset of Fig- El we plot 
the dependence of S y (t) / S y (0) against t for two cases by 
the solid and dashed curves. In these two cases 5^(0) 
are initialized from the intrinsic situation at 150 K under 
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FIG. 8: (Color online) Spin relaxation time r 3 against the 
initial spin polarization S y (0) for the intrinsic (solid curves) 
and n-type (dashed curves) cases at different temperatures. 
The temporal evolution of S y (t) / S y (0) is plotted in the inset 
for two cases, with S y (0) initialized from the intrinsic situa- 
tion at 150 K under electric field 1 (solid curve) and 7 kV/cm 
(dashed curve), respectively. The dotted (chain) curve is the 
exponential fit to the rapid decay in the beginning of the solid 
(dashed) curve. 



electric field 1 (solid curve) and 7 kV/cm (dashed curve), 
respectively. To obtain the spin relaxation time r s , we 
fit the rapid decaying part in the beginning by an expo- 
nential function, as illustrated by the dotted and chain 
curves in the inset. 

The spin relaxation time r s against the initial spin po- 
larization S y (0) for the intrinsic and n-type cases at dif- 
ferent temperatures is plotted in Fig. [8j It is shown that 
the spin relaxation time is of the order of 0.01-0.1 ps, 
mainly limited by the electron-phonon scattering. In 
fact, a momentum relaxation time of the similar order 
is given from both estimation and experiment by Butch 
et al.M- The decrease of r s with the increase of S y (0) 
is due to the fact that the larger S y (0) is initialized by 
the higher electric field, with which the electrons and 
holes occupy larger-momentum states and hence exhibit 
faster momentum relaxation. Besides, with higher tem- 
perature, the spin relaxation time also decreases as the 
electron-phonon scattering is strengthened. 



V. CONCLUSION 

In conclusion, we have studied the charge and spin 
transport under the influence of high electric field (up 
to several kV/cm) on the surface of topological insula- 
tor Bi2Se3, by means of the KSBEs. We assume that 
the Fermi level, adjustable by doping^&22r— is located 
in the bulk gap. Therefore the bulk states are excluded 
from our study. With moderate electron and hole densi- 
ties, the surface state around the Dirac point can be de- 
picted by the Rashba spin-orbit coupling. In our study, 



both the conduction and valence bands of the surface 
state are considered, with the inter-band coherence ex- 
plicitly included. Apart from the driving effect in each 
band, due to the spin mixing of conduction and valence 
bands, the electric field also leads to inter-band preces- 
sion. This differs from the semiconductors with parabolic 
energy spectrum. 

In Bi2Se3, the relative static dielectric constant is 
as large as 100, indicating weak electron-impurity and 
electron-electron Coulomb scatterings. With the weak 
Coulomb scattering, electrons in the two bands fail to 
establish a drifted Fermi distribution with a unified hot- 
electron temperature under the driving of the electric 
field. The electron-surface optical phonon scattering 
dominates in a large temperature region. Moreover, the 
electron-phonon scattering matrix element is approxi- 
mately constant due to its marginal dependence on mo- 
mentum. This feature leads to particular properties of 
charge and spin transport on the surface of Bi2Se3. 

Our study reveals that in the presence of driving of the 
electric field, both the inter-band precession and inter- 
band electron-phonon scattering cause electrons to trans- 
fer from the valence band to the conduction one. Due to 
the dominant electron-phonon scattering, the variation 
in electron density for each band is linear in the elec- 
tric field when the latter is high, despite whether the 
initial state is degenerate or nondegenerate. It is also 
found that due to the spin-momentum locking from the 
Rashba spin-orbit coupling, a transverse spin polariza- 
tion is induced by the electric field, with the magnitude 
proportional to the momentum scattering time. Besides, 
the spin polarization is linear in the electric field when 
the latter is small but deviates from the linear relation 
when the latter is large enough as the electron-phonon 
scattering is enhanced due to the heating of electrons. 
Moreover, a very interesting feature is that at high tem- 
perature, the spin polarization is inversely proportional 
to the temperature but insensitive to the electron density. 

The cooling of hot carriers and the relaxation of spin 
polarization induced by the electric field are investigated 
by turning off the electric field after reaching the steady 
state. It is found that the hot carriers cool down in a 
time scale of energy relaxation, which is quite long and 
is of the order of 100-1000 ps. However, due to the spin- 
momentum locking again, the spin polarization relaxes 
in a time scale of momentum scattering. The spin po- 
larization is mainly contributed by the states with large 
momentum. Therefore, it decays rapidly within the time 
of the order of 0.01-0.1 ps. Following this rapid decay, 
there is a slowly damping tail. This tail is attributed 
to the low-energy states near the Dirac point where the 
momentum scattering is weak. 
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withS k (>,<) = S?(>,<) + Sg(>,<) + S?{>,<) from 
the electron-surface optical phonon, electron-impurity 
and electron-electron scatterings, respectively. Here 



Appendix A: Scattering term in KSBEs 

The scattering term in Eq. @ can be written as 

dtPkWIscat =S k (>, <) - S k (<, >) + St(>, <) 

-St(<,>), (Al) 



^k P (>, <) =| \9^ 2 ^M'Pw(t)T m A k >MT V2 p<(t)[N<S{e klm - e kl?2 - hw ) + N>S(e^ Vl - e kl?2 + (A2) 
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(A3) 



(A4) 



To obtain the above equations, the Markovian approx- 
imation, Pk(t') ~ p\J\t)i is adopted in the interaction 
picture. Here p\Jt) = e iEkt pk(t)e~ kt with the matrix 
Ek = diag{e k +, e k _}. Equivalently, in the Schrodinger 
picture, the Markovian approximation is written as 
Pk(t') « e-*®^' -V pk{t)e iB *& -V . With this approxima- 
tion, the scattering terms of KSBEs2£ are consistent with 
those of kinetic Bloch equations in semiconductors^ 
However, in the work by Culcer et al, the Markovian ap- 
proximation is mistakenly carried out as Pk(t') « Pk(t) 
in the Schrodinger pictured Therefore, the scattering 
term given there deviates from ours. In Eqs. (|A2[) - (|A4|1 . 



A 



k,k' = u£u k > [U k is given in Eq. ©] and T m = ±(1 + 
jjiffg) with r)i = ±. p£(t) = pk(t) and p£(t) = 1 - Pk(t). 
Nq = N + 1 and AT< = N where N is the Boson dis- 
tribution of surface optical phonons with energy hioo = 
7.4 meV.— |g q | 2 is the electron-phonon scattering matrix 



[Uk is given in Eq. ©] and T v 



\D^ = \ZiV q °e 



- 9rf /e(q,0)| 2 is 



element given by Eq. ([7} . 

the electron-impurity scattering matrix element, where 
Zi is the charge number of impurity (assumed to be 1 
here), d is the effective distance of impurities from the 
surface two-dimensional electron layer (assumed to be 
zero following Culcer et alM-), and V® = 2nhvfr s /q with 

r s = e 2 /(4 Wo ^/). In Eq. flU, V£(fej) = V?/e(q,w) 
and Vq 2 (hcu) — Vq*(huj) are the retarded and advanced 
Coulomb potentials with dynamic screening, respectively. 



The RPA screening e(q, w) = 1 - V g n(q, with 



n(q,w) = V t /km /k+q " 2 

— ' nuj + £ kj)1 — e k+ql)2 
x [1 + 771/72 cos(6» k - 6» k+q )]/2 



i0+ 



(A5) 



In the n-type degenerate (intrinsic nondegenerate) sit- 
uation, under the long-wavelength and static limit, 
11(0,0) = -kf/(2irtiVf) (Ref. M) [n(0,0) = 
-k B T\nA/(2TTh 2 vj)} given by Eq. (|A"5t with / kl; substi- 
tuted by the equilibrium Fermi distribution. Therefore 
e(0,0) = 1 + q s /q with q s = r s k f (q s = r s k B TlnA/hvf) 
for the n-type degenerate (intrinsic nondegenerate) case. 
Especially, for the intrinsic case, due to the large Fermi 
velocity Vf and small r s , the screening is very weak and 
we further approximate e(0,0) ~ 1. These approxima- 
tions for screening are utilized in analytically calculating 
the momentum relaxation time limited by the electron- 
impurity scattering in Sec. IIII Al 



Appendix B: Energy loss rate due to the 
electron-phonon scattering 

The energy loss rate, induced by the electron-phonon 
scattering solely, satisfies 



d t s (t) = J2 hv f kTr[a z d t Pk(t)\2j. (Bl) 
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Near the equilibrium, approximating pk(t) isotropically 
and retaining its diagonal part only, one has 



d t e (t) 



4nh 



■jNfco/ dkk(k + k )[N>f% +ko 



Jt) 



fro 



- N Q <r k (t) - r k {t)r k+ko {t)\ + i dkk\k v 

Jo 

x + # (*) - 1) + ns)frm}- 

(B2) 

Here = / k+ (i) and /^(t) = 1 - / k _(*) stand for 

electron and hole distributions in the conduction and va- 
lence bands, respectively. On the right-hand side of the 
equation, the first integral is contributed by the intra- 
band electron-phonon scattering in both bands while the 
second one by the inter-band electron-phonon scattering. 
For the n-type degenerate case with both low tempera- 
ture and electron density, we take into account the con- 
tribution of the intra-conduction band electron-phonon 
scattering to calculate the energy loss rate as 



dt£o(t) 



M 2 fcp 2 



dkk[N>f+ +ko (t)-N<f+(t) 



-fk(t)fk +k0 (t)} 



(B3) 



By using the detailed balance condition satisfied in 
the initial equilibrium state, 7V > f£ +ko (0) - N^f£(0) - 

/^(O)ZaWO) = 0, one has 



9te (t) £°° dkk{N>[f+ +ka (t) f+ +ko (0)] 

-Kitiv-rtm-fitwz+ko® 

+ / fe + (o)/ fc + +feo (°)} 



lffo| 2 fcg 

Anh J 

\9o\ 2 k* 



p-\-oo 



dkk\f+(t) - /+(0)] 



2h 



n e (t)-n e (0)]. 



(B4) 



This relation is utilized to obtain Eq. ([T 
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